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INTRODUCTION : 


/ 

"  When  a  body,  floating  on  the  surface  of  an  infinite, 
ideal,  inviscid,  irrotational  fluid  is  subjected  to  a  periodic 
vertical  displacement,  a  wave  pattern  is  created  in  the 
fluid  and  the  problem  of  determining  this  pattern  from  a 
knowledge  of  the  body  geometry  and  applied  forces  is  well 
known  in  fluid  mechanics.  In  problems  with  both  partially 
and  fully  submerged  objects,  quantities  of  physical 
interest  are  not  only  the  wave  patterns  which  may  be  derived 
from  the  velocity  potential  but  also  functionals  of  the 
potential  such  as  added  mass  and  damping  factors  which 
measure  the  distribution  of  energy  in  the  fluid,  e.g. 

Weyhausen  and  Laitone  [26,  p.  567].  These  factors  are,  of 
course,  dependent  on  the  body  geometry . and  the  present  paper 
is  devoted  to  showing  how  these  quantities  may  be  optimized 
over  restricted  classes  of  body  geometry. 

Specifically  we  Will  study  the  problem  of  the  optimal 
design  of  a  floating  body,  totally  submerged  in  a  fluid  of 
finite  depth.  In  the  terminology  of  optimal  control,  this 
is  a  problem  of  optimization  of  geometrical  elements 
(see  Lions  [17 ] ) . 

In  his  classic  paper  [13],  John  formulated  this  problem 
of  a  partially  immersed  heaving  body  as  a  boundary  value 
problem  for  the  velocity  potential  which  satisfies  Laplace's 
equation  with  given  Neumann  data  on  the  submerged  portion 
of  the  boundary,  a  linearized  free  surface  condition  on  the 
mean  free  surface,  a  homogeneous  Neumann  condition  on  the 
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bottom  of  the  fluid  container,  and  a  radiation  condition. 
Under  certain  restrictions  on  the  body  shape,  in  particular 
that  it  be  convex,  smooth,  and  have  normal  intersection  with 
the  free  surface,  he  proved  the  existence  of  a  unique 
solution  and  formulated  an  integral  equation  for  the 
velocity  potential  evaluated  on  the  body  surface.  Kleinman 
has  shown  in  [15]  that  the  geometrical  conditions  of  John 
may  be  somewhat  relaxed:  corners  are  allowable  as  are 
non-normal  intersections  with  the  free  surface,  and  convexity 
is  not  necessary.  However  the  essential  restriction  that 
vertical  rays  from  the  free  surface  intersect  the  body  at 
most  once  remains.  Relaxation  of  this  restriction  has  been 
reported  by  Ur sell  [25]  for  the  two  dimensional  case.  In 
[15]  it  is  also  shown  that  the  integral  equation,  suitably 
modified  possesses  unique  solutions  i.e.  methods  are 
provided  for  showing  how  one  may  eliminate  the  problem  of 
the  existence  of  so-called  irregular  frequencies.  We 
should  point  out  here  that  such  irregular  frequencies  are 
an  artifact  of  the  integral  equation  method;  they  are 
particular  values  of  the  parameter  appearing  in  the  free 
surface  boundary  condition  for  which  the  integral  equation 
has  multiple  solutions.  However,  the  original  boundary 
value  problem  has  at  most  one  solution  and  so  we  see  that 
an  irregular  frequency  does  not  correspond  to  a  physical 
pathology.  Moreover,  such  irregular  frequencies  are 
dependent  on  the  shape  of  the  body.  As  we  will  be  looking 
at  a  family  of  such  bodies  it  will  be  essential  to  our 


analysis  that  we  use  an  integral  equation  formulation  that 
is  free  of  such  irregular  frequencies  and  one  of  our  first 
tasks  will  be  to  establish  such  a  boundary  integral 
equation  for  our  problem. 

When  the  body  is  completely  submerged,  John's 
uniqueness  proof  no  longer  applies,  however,  Maz'ja  [19] 
has  provided  a  proof  for  a  class  of  bodies  subject  again 
to  certain  geometric  restrictions.  In  this  case  the 
boundary  value  problem  may  once  more  be  reformulated  as  a 
uniquely  solvable  integral  equation.  A  virtue  of  the 
integral  equation  formulation  of  problems  for  partial 
differential  equations  in  unbounded  regions  exterior  to 
a  bounded  obstacle  lies  in  the  apparent  advantage  that 
numerical  algorithms  for  bounded  domains  (the  domains  of 
the  integral  operators)  have  over  algorithms  for  unbounded 
domains. 

Problems  involving  optimization  with  respect  to  the 
domain  have  been  treated  by  a  number  of  authors,  among 
whom  we  mention  Cea  and  his  co-workers  [  5  ] ,  [  6  ] , 

Chenais  [7  ]  and  Pironneau  [21],  [23].  This  last  author 

has  published  a  book  [24]  on  this  general  topic.  But  the 
great  majority  of  the  papers  in  this  area  treat  problems 
in  bounded  domains.  In  this  paper,  we  deal  with  a  problem 
whose  natural  setting  is  an  unbounded  domain  and  we  will 
find  the  boundary  integral  equation  method  particularly 
convenient.  With  respect  to  optimization  problems  in 
exterior  domains,  using  boundary  integral  equations,  we 


refer  to  the  recent  treatment  of  the  three  dimensional 
inverse  acoustic  scattering  problem  by  Angell,  Colton,  and 
Kirsch  [  1  ]  and  the  subsequent  numerical  treatment  of  an 
inverse  two  dimensional  problem  by  Kirsch  [14] . 

As  in  [  1  ] ,  our  approach  to  the  optimization  problem 
will  depend  on  the  reformulation  of  the  original  boundary 
value  problem  (which  here  includes  not  only  boundary 
conditions  given  on  the  bounded  surface  of  the  body,  but 
also  on  the  bottom  and  on  the  free  surface  which  are  of 
infinite  extent) ,  as  a  uniquely  solvable  integral  equation 
defined  on  the  boundary  of  the  body.  The  first  section  of 
this  paper  is  devoted  to  this  question.  While  the  use  of 
an  appropriate  radiation  condition  and  Green's  theorem 
will  lead  to  an  integral  equation,  it  is  important  to  note 
that  the  equivalence  of  this  integral  equation  formulation 
and  the  original  boundary  value  problem  depends  on  a 
uniqueness  theorem  for  the  boundary  value  problem.  This  is 
already  a  non-trivial  problem  which  is,  as  yet,  not  fully 
understood.  As  mentioned  earlier,  for  the  case  of  a 
totally  submerged  body,  Maz'ja  [19]  has  given  a  proof  of 
uniqueness  for  a  restricted  class  of  bodies.  The  recent 
and  interesting  paper  of  A.  Hulme  [21]  discusses  the 
result  of  Maz'ja  and  most  effectively  describes  the 
geometric  meaning  of  the  result.  Since  it  is  essential  to 
have  a  uniqueness  result  for  the  original  boundary  value 
problem,  we  will  confine  ourselves  to  this  class  of  bodies; 
our  results  can,  however,  be  extended  immediately  to  any 
situation  in  which  uniqueness  obtains. 


The  second  section  of  the  paper  is  then  devoted  to 
the  formulation  and  proof  of  the  existence  of  a  solution, 
optimal  in  the  sense  of  being  a  body  with  minimal  (or 
maximal)  added  mass  or  damping,  in  any  compact  set  of 
admissible  surfaces.  We  point  out  here  that  this  problem 
is  highly  non-linear:  the  cost  function  involves  the 
solution  of  the  boundary  value  problem  relative  to  the 
unknown  boundary. 

Since  the  class  of  all  admissible  surfaces  is  convex 
we  study  necessary  optimality  conditions  in  section  3. 

There,  we  compute  the  Gateaux  derivative  of  the  cost 
functional  using,  basically,  the  arguments  introduced  by 
Hadamard  [11],  The  problem  of  computing  the  derivative  of 
a  "domain  functional"  appears  in  a  number  of  works 
including  [2],  [5],  [6],  and  [  7  ]  .  The  techniques  used 

here  allow  us  to  deal  with  a  much  more  general  situation 
than  that,  for  example,  of  Kirsch  [14]  in  so  far  as  we 
work  in  three  dimensions  and  make  no  a  priori  analyticity 
assumptions  on  the  parametric  representation  of  the  surface. 
Moreover  we  do  not  use  techniques  of  complex  variables  e.g. 
conformal  mapping  arguments  or  reflection  principles. 

As  one  would  expect,  since  we  are  dealing  with  a 
constrained  optimization  problem,  the  necessary  conditions 
take  the  form  of  a  variational  inequality  and  we  conclude 
section  3  with  a  brief  discussion  of  a  numerical  procedure 
for  finding  solutions  of  this  inequality.  It  is  our 
intention  to  present  concrete  numerical  results  elsewhere. 
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1 .  The  Basic  Boundary  Value  Problem 

In  this  discussion  we  will  consider  the  following 

geometry:  A  Cartesian  coordinate  system  is  centered  at  a 

point  on  the  linearized  free  surface  and  oriented  so  that 

the  free  surface  coincides  with  the  (x,z) -plane  (y=0)  and 

the  fluid  extends  from  y=-h  to  y=0 .  The  submerged  body 

will  be  a  simply  connected  subset  of  R  lying  in  a  strip 
2 

R  xf-h,-^]  whose  boundary  r  is  a  bounded  Lyapunov 
surface  of  index  1.  We  will  place  certain  additional 
restrictions  on  F  and  these  will  be  indicated  as  they 

become  necessary.  We  denote  the  interior  of  the  body  by 

—  2 

D  ,  and  the  points  of  the  strip  R  x[-h,0]  and  exterior  to 

the  body  by  D+.  The  condition  that  the  surface  be 
Lyapunov  of  index  1  guarantees,  among  other  things,  that 
there  exists  a  Lipschitz  continuous  normal  n  at  all  points 
of  r.  We  emphasize  that  n  is  oriented  so  that  it  points 
into  D+.  Points  will  be  denoted  by  p=(x  ,y  , z  )  with 

r  r  r 

cylindrical  coordinates  p=(p  , 0  , y  )  and  the  subscripts 

P  P  P 

will  be  omitted  if  there  is  no  danger  of  confusion. 


We  concern  ourselves  with  the  following  boundary  value 


problem: 


(1.1) 


(a) 

o 

-©- 

n 

o 

in  D 

(b) 

In  + 

=0  on  y=0 

(c) 

0 

on  y=-h 

.(d) 

on  r 

together  with  a  radiation  condition 

(1.1)  (e)  |i  -  ik04>  =  o{p~H) 

where  gee (F)  and  kQ  is  the  root  with  largest  real  part  of 
the  transcendental  equation 


(1.2)  k  sinh  k  h  *  k  cosh  k  h  . 
n  n  n 

Recall  that  n  always  points  into  D+ . 

In  [19],  Maz'ja  introduced  a  restricted  class  of 
boundaries  for  which  this  boundary  value  problem  has  at  most 
one  solution.  We  formulate  that  theorem  here  as  follows: 


Theorem  1.1:  Let  V  be  the  vector  field  in  defined  by 


V 


P( 


2  2, 
-P  ) 


p  +y 


P  - 

r  2  i 

2p  y 

2,  2 

J 

Lp  +y  . 

Then  the  homogeneous  boundary  value  problem  (1)  with  g=0 
has  only  the  trivial  solution  provided 


(1.3)  V*n  >0  on  F 
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A  discussion  of  this  result  and  its  geometric  significance 
may  be  found  in  A.  Hulme  [12] .  We  will  refer  to  the  class 
of  all  such  surfaces  as  the  Maz'ja  class. 

Following  John  [13]  we  introduce  the  Green's  function 
for  this  problem  which  is  normalized  to  have  the  form 

(1.4)  y(p,q)  =  +  R(p,q) 

where  the  function  R  has  bounded  derivatives  with  respect 
to  q  for  points  per  (see  [13;  p.  96])  and  y  satisfies 
conditions  1.1b,  c  and  e.  Using  this  Green's  function  to 
define  single  and  double  layer  potentials,  the  usual  jump 
conditions  can  be  established  as  in  the  potential-theoretic 
case  since  the  singular  behavior  of  y  and  •£— -y  is  determined  by 

■  nq 

the  first  term  in  (1.4) .  For  convenience,  we  record  these 
results  here; 


(1.5)  lim  u(q)y(p,q)dr  =  ±u(p)+  u(p) 

o+r-  dnn  3rV, 


dr  , 
q 


(1.6)  lim+  [  u  (q)  ~  Y(p,q)dr  =  +u(p)  +  f  u(q)  aYjP,q)  dr 

p->r_  j,  °q  q  f  dnq  q 

where  p-*T_  means  p  approaches  T  from  D*. 

Moreover,  if  u  is  a  solution  of  the  boundary  value 
problem  (1.1)  then  one  may  use  Green's  Theorem  to  draw  the 
familiar  relation 

,  r  1  C2u(p),  peD+ 

(1.7)  j  Y(p,q)  (y(p,q,  )u(q)  |dr  =  ■{  u(p)  ,  per 
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If  one  then  uses  the  boundary  condition  (l.ld)  we  have, 
(1.8)  Y(p,q)g(q)drg-j  u  (q)  [y(p,q)]drq  =  u(p) 


or,  in  operator  notation 


(1.9)  (I+K*) 


f 

u  = 


y (p,q) g (q) dr 


where  K*  is  the  boundary  integral  operator  with  kernel 

q 

We  pause  to  remark  that,  given  a  solution,  w,  of  this 
integral  equation  we  may  represent  the -solution  of  the 
boundary  value  problem,  according  to  the  relation  (1.7)  by 

(1.10)  u  ( p)  =  hi  y  (p,q)g(q)dr  -%f  w(q)  dr  ,  peD+ 

q  j  on  q 


and,  again  using  the  jump  relations  one  sees  easily  that 


(1.11)  u  r  =  u>  , 

which  is  a  direct  relationship  between  the  solution  of  the 
boundary  integral  equation  and  the  boundary  values  taken  on 
by  the  solution.  Such  a  direct  relation  does  not  obtain 
when  one  uses  a  layer  approach  in  which  one  assumes  that 
the  solution  u  has  a  representation  as  a  single  layer. 


u ( p)  =  u>(q)Y  (p,q)dr  , 


and  then  the  boundary  condition  and  jump  relations  are 
employed  to  obtain  an  integral  equation  for  co. 

Our  first  task  is  to  discuss  the  unique  solvability 
of  (1.9) . 


Theorem  1.2:  Let  r  be  Lyapunov  of  index  1  and  belong  to 
the  Maz'ja  class.  Let  geC(r).  Then  the  integral  equation 
(1.9)  has  a  unique  solution. 

Proof :  Since  K*  is  an  integral  operator  with  a  weak 

singularity  on  the  smooth  surface  r  it  is  compact 

(e.g.  [  8  ])  hence  the  Fredholm  alternative  applies.  Thus 

the  integral  equation  in  question  either  has  a  unique 
2 

solution  in  L  ( r  )  (or  in  fact  in  C ( T  ) )  or  the  homogeneous 

adjoint  equation  has  non-trivial  solutions  in  which  case 

f 

(1.9)  has  solutions  if  and  only  if  the  function  ygdF  is 

orthogonal  to  all  solutions  of  the  homogeneous  adjoint 
equation 

(1.12)  (I+K)S  =  0  , 

where  K  is  the  boundary  integral  operator  with  kernel  . 

np 

We  first  prove  that  any  solution  to  (1.12)  is  orthogonal 
to  |  ygdr  for  any  geL^(T). 

r 

£ 

Assume  that  u  is  a  solution  of  (1.12)  and  define 

r  .  £  + 

<MP)  :=  y  (p,q)  u)  (q)  dF  ,  peD 

£ 

Certainly  w,  as  a  solution  of  (1.12)  is  continuous  on  F 
and  the  jump  condition  (1.5)  yields 

=  w+Koo  =  0  on  F  . 

Hence  the  function  \p  is  a  solution  of  the  exterior 
homogeneous  submerged  floating  body  problem  and  so,  by 


Theorem  1.1,  ^=0  in  D' .  Therefore, 

0  =  |  g  (q)  (q)  dig  =  j  Y  (q,p)  <2  (p)  drpg  (q)  drq 

r  r  r 

and  because  the  Green's  function  is  symmetric  we  obtain 

j  Y(p,q)g(q)drqiS(p)drp  =  o 
r  r 

which  is  the  required  orthogonality  condition.  Thus  the 
equation  (1.9)  always  has  at  least  one  solution. 

To  see  that  (1.9)  has  at  most  one  solution,  we  assume 

O  *  _  £ 

that  w*  is  a  non-trivial  solution  of  (I+K*) S*=0 .  Define  v 
in  D_  by 

v(p)  =  f  LY<j?'Sl  s*(q>ar  ,  P£D‘  . 

Then,  again  using  the  jump  relations  we  have 

v(p)  =  (I+K*)C*(p)  =  0 

and  so  v  is  a  solution  of  the  homogeneous  interior  Dirichlet 
problem  for  the  Laplacian.  Since  there  are  no  non-trivial 
solutions  of  this  problem,  v  must  vanish  in  D“  and  so 
( p)  =0  for  per 

9n~ 

Now  define  v  in  D+  by 

v(p)  =.  f  ^P-3.1-  5*(q)ar  ,  P£D+  . 

]  3nq  q 

o 

Since  Jj*  is  continuous  on  r  and  r  is  Lyapunov  of  index  1, 

~~ (p )  exists  and  is  equal  to  -^-(p)=0  by  a  Theorem  of 
3n+  9n“ 

Lyapunov  (Gunter  [  9 ?p.  297]). 


So  again  by  Maz'ja's 


uniqueness  theorem,  v  must  vanish  in  D  and  so,  by  the 
jump  conditions  we  have 

-u*+K*5*  =  0  . 

O  _  O  o 

But  S3*+K5*=0  as  well,  hence  S*=Q  which  is  a  contradiction 
and  the  theorem  is  established. 


L3 


2 .  The  Optimization  Problem 

*3  I 

Let  fg={peE  |p|=l}  denote  the  surface  of  the  unit 
ball  in  E3  and  let  C1'^(ro)  denote  the  space  of  continuously 
differentiable  functions  whose  first  derivatives  satisfy  a 
Lipschitz  condition  and  which  is  equipped  with  the  usual 
Holder  norm  | )• | |^  ^  (see  e.g.  [  8  ]) .  We  will  assume  that 
we  are  given  a  family  of  surfaces  which  can  be  described 
by  parameterizations : 


(2.1)  r (f)  =  {p£EJ | P=f (p) p+pg,  p  =  |p—  --j} 

where  f:rQ->-E3  is  an  element  of  and  pQeE2x  (-h,-€Q)  . 

Let  a  and  b  be  two  positive  constants  and  define  the  subset 

Fa,bccl,l(ro> 

(2*2)  Fa,b  =  {fecl'1(r0)  I  I  if  I  llfl<b,  f (p)p+P0eE2*(-h,-e0) 

and  f (p) £a,  PerQ} 

.  .  3 

Definition  2.1;  A  surface  S  in  E  will  be  called  admissible 
provided  S  can  be  described  by  a  parameterization  f«F=  ,  and 
S  is  contained  in  the  Maz' ja  class  (c.f.  Theorem  1.1) . 

Note  that  since  each  admissible  surface  is  completely 
determined  by  the  function  f,  we  will  henceforth  simply 
refer  to  "the  surface  f"  although,  when  convenient,  we  will 
use  the  notation  T(f).  Clearly,  each  admissible  surface 
describes  a  surface  bounding  a  bounded  region  which  contains 
a  ball  of  radius  a/2  and  center  pft  in  its  interior. 
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Figure  2 

,2. 


We  will,  when  necessary,  denote  the  region  in  £  *(-h,0) 
exterior  to  an  admissible  surface  f  by  D+  and  that  interior 
to  the  surface  by  . 

Let  Uad  be  a  compact  subset  of  admissible  surfaces. 


2  11 

For  example,  since  the  imbedding  C  (rQ)-+C  '  (rQ)  is 


compact,  we  may  choose  Uad  to  be  the  subset  of  admissible 
,  2 

functions  in  C  ( r Q )  .  This  particular  choice  leads  to  a 
non-linear  optimization  problem  over  a  closed  convex  set. 


The  convexity  will  be  advantageous  for  subsequent  numerical 
considerations . 

We  now  wish  to  consider  a  family  of  boundary  value 
problems  of  the  type  discussed  in  section  1  which  may  be 
considered  as  indexed  by  uad: 


(2.3) 


ri  a) 

A <J>  (p) 

= 

0,  peD*  , 

(b) 

3(j> 

Jn  + 

k<p 

=0  on  y=0 

{  (c) 

3<J> 

TE  ~ 

0 

on  y=-h  , 

1  <d) 

d<p 

7n  ~ 

g 

on  T (f )  , 

L(e) 

3  <£ 

Tp  " 

ik 

_  it 

0<p  =  o(p  . 

Note  that,  because  we  are  considering  a  family  of  boundary 
value  problems,  the  data  g  in  (2.3d)  must  be  defined 
throughout  the  domain  formed  by  the  union  of  all  admissible 
surfaces.  This  is  indeed  the  case  for  heaving  motion 
where  g=-n*y. 

With  this  understanding,  each  choice  of  surface 
f£Uad  gives  rise,  according  to  Theorem  1.1,  to  a  potential 
<p=<p(p;f),  p^Dj.  The  class  of  optimization  problems  that 
we  discuss  below  then  have  the  form:  Let  J:C{F0)->-E  be 
continuous  and  define  a  map  J  [  •  ]  : U  .-►!  by 

(2.4)  J[f]  :=J($  ( - ;  f )  )  ,  feUa{J 
where 

(2.5)  $  ( p?  f)  :=4>  (f  (p)  P+P0;  f )  ,  P*r0  . 

The  optimization  problem  is  then  to  find  f0eUad  suc^  that 

(2.6)  J[fQ]  <  J(f],  for  all  feUad 
or 

(2.7)  J[fQ]  >  J[fJ,  for  all  feUad  . 

We  will  confine  our  discussion  to  the  problem  of  minimization 
This  is  sufficiently  general  since  the  problem  of  maximizing 
a  functional  J  may  always  be  replaced  by  that  of  minimizing 
-J. 

Specific  forms  of  the  functional  J  of  (2.4)  may  be 
chosen  to  reflect  desirable  design  criteria.  For  example, 
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as  mentioned  in  the  introduction  one  may  choose  J  to 
represent  the  added  mass  of  the  hull.  In  this  case,  the 
problem  of  interest  is  that  of  minimizina  the  functional 
J  in  order  to  reduce  the  hydrodynamic  force  on  the  ship 
hull,  a  goal  which  is  obviously  of  great  importance  to 
ship  design.  Indeed,  it  is  well  known  (see  [26;pp.  563-567]) 
that  the  added  mass  of  a  particular  hull  may  be  represented  by 

Ma  *  Re  l  *<P>  &  • 

r  (f ) 

This  form,  in  light  of  the  boundary  condition  (2.3d),  leads 
to  the  functional 

* 

(2.8)  J[f]  =  Re  $(p;f)g(f (p)p+pQ)Jf (p)dr~  , 

'  i 

0 

where  Jf  is  the  Jacobian  of  the  transformation  p=f(p)p+pQ. 

As  we  will  see  below,  this  form  indeed  defines  a  continuous 

functional  on  U  , . 

ad 

Returning  to  the  more  general  problem  of  minimizing 
the  functional  (2.4),  our  first  task  is  to  show  that  the 
functions  $(*;f)  related  to  the  family  of  boundary  value 
problems  (2.3)  via  the  equation  (2.5)  depend  continuously 
on  the  choice  of  surface  f,  i.e.  that  the  mapping 
f**$  ( • ;  f )  is  continuous  as  a  map  from  Uad  into  C(r0)  . 

Recall  that  the  solution  of  the  boundary  value  problem  (2.3) 
is  given  in  (see  also  equation  (1.10))  by 

(2.9)  <p  (p)  =  Jjf  Y(p,q)g(q)dT  -hf  1—  Y(p,q)u(q)dl\  , 

q  dna  q 

r (f )  r (f )  q 


where  u  is  a  solution  of  the  integral  equation  (1.9)  on 


r  ( f )  .  Thus  <Hp;f)|  =u(p;f)  and  we  may  write  the 

1  r  (f) 

integral  equation  (1.9),  suppressing  the  dependence  of  <J> 
on  f  as 


(2.10)  <f 


y  (p,q)  <p  (q)drq 


Y(p,q)g(q)drg  . 
(f) 


Theorem  1.2  guarantees  that  this  equation  possesses  one 
and  only  one  solution  for  each  given  feUad  and  geC(T(f)). 

We  now  wish  to  use  the  integral  equation  (2.10)  to 
study  the  mapping  f^$(«;f)  as  a  map  from  Uad  into  C(rQ) . 

In  particular  we  wish  to  show  that  this  map  is  continuous 
and,  consequently,  a  functional  J  of  the  form  given  by  the 


relation  (2.4)  will  assume  its  minimum  value  on  U 


ad* 


To 


this  end,  we  consider  the  integral  equation 


(2.11)  d>(f  (P)P+P0)  + 


•  3 


[•jpf"  Y(f  {p)p+pQ,f  (q)q+pQ)  ]<Mf  (q)q+pQ)  Jf  (q)dr 

q 


=  f 


j  Y (f (p) P+P0rf (q)q+P0)g(f (q)q+P0) Jf (q)df~ 


obtained  from  the  integral  equation  (2.10)  by  means  of  the 
change  of  variables  p=f(p)p+pQ.  Note  that,  under  our 
hypotheses,  the  Jacobians,  J^,  of  the  transformations,  are 
uniformly  bounded  on  rQ  with  respect  to  f£Ugd  since  Uad  is 
a  bounded  set  with  respect  to  the  Clfl-norm.  Indeed,  writing 

q=f  (6q,4>q)q+P0,  q=(sin6q  cos<f>q,  sinGq  sin<J)q,  cosGq)  ,  then, 

for  any  given  f£Ua(j»  the  Jacobian  has  the  form 


(2.12)  J*(6„,4>  )=f  (9  ,6  )  [f  ( 9_ ,  <t>  )  sin9  +(^f(6  .4  )) 


f v  q'yq 


q'  rq,  «•*  '  q  q 


&ia  °q  M  a'V' 

+  (JF  f  Oq-tg)  >  2sin2e<3]ls. 


Since  each  term  of  the  second  factor  on  the  right  is  bounded 
by  I  I  f  I  1 1,  i  and  lf  (eq/<J>q)  I- I  !f  I  1 1,1  we  see  that 

(2.13)  sup  | Jf (q)  j  £  /?  ( | f } [ ^  1  Z  b2/l  . 

q^r0 

Introducing  the  functions  $  (p)  :=<j>  (f  (p)  p+pQ)  * 
gf  (p)=g(f  (P)P+P0)  ,  af  (p,q)  =  [|jj-  Y(f  (p)p+p0#  f  (q)q+P0)  1  jf  (q) 

q 

and  bf (p,q)=[y(f (p)p+PQ;f (q)q+PQ) ] Jf (q) #  then  equation  (2.11) 
assumes  the  convenient  form 

(2.14)  $(p)+  af (p,q)$(q)dr~  =  bf (p,q) gf (q) df-  . 


Let  B(C(r0) )  denote  the  class  of  bounded  linear  operators 
on  C (Tq) .  We  define  two  families  of  operators  (Af}feU 

ad 

and  in  8(C(rQ))  by  the  relations 

ad 

* 

(2.15)  Af  (ip)  (p)  =  af  (p,q)  4»(q)  dT~,  p*rQ 


(2.16)  Bf(*)(p)  =  bf  (p,q)U>(q)dr5,  P^( 


We  will  establish  the  continuity  of  the  map  fH-$  by  showing 
that  the  kernels  a^  and  b^  are  weakly  singular  for  all 
f<Uad  and  that  the  operators  A^  and  depend  continuously 
on  f  i.e.  that  the  maps  fn-A*  and  fH*-B*  of  Ua^  into  8(C(r„)) 
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are  continuous  when  8(C(rg))  is  equipped  with  the  uniform 
operator  topology.  To  this  end,  we  will  need  two  crucial 
estimates  which  are  given  in  the  following: 


Lemma  2.1:  For  any  6,  0<6<Js,  the  kernels  af  and  bf  are 
weakly  singular  and  satisfy  the  inequalities: 

c, 


(2.17)  | af (p, q) -a  (p, q) |  <  | | f-g 


<5 


1,1  |p-q 


T+7S 


and 

(2.18)  |b  (p,q)-b  (p,q)  |  <  |  |  f-g  j  |  6  - Kt-t 

9  1,1  |p-q|1+6 

where  c^  and  c2  are  positive  constants.  The  proof  is 

quite  technical  and  we  relegate  it  to  an  appendix.  An 

estimate  similar  to  (2.17)  was  used  in  [  1  ].  The  present 

estimates  are  used  in  the  proof  of  the  first  part  of  the 

next  theorem. 


Theorem  2.1:  Let  8(C(Fg))  denote  the  space  of  bounded 
linear  operators  on  C(Tq)  equipped  with  the  uniform 
operator  topology  and  assume  that  the  map  fn-g^.  from 
C1'1^)  into  C(Tq)  is  Holder  continuous.  Then  the  mappings 
fKAf  and  f*^Bf  of  C1,:L(r0)  into  B(C(Tq))  are  Holder  continuous. 
Moreover,  since  the  set  Uad  is  compact,  the  map  fn-$(.;f) 
of  Uad  into  C ( rQ) ,  where  $(*;f)  is  the  unique  solution  of 
the  boundary  integral  equation  (2.14),  is  Holder  continuous. 

Proof :  Consider,  first,  the  mapping  f^Bf.  Using  the 
estimate  (2.18)  we  have 


ffl 


Bf  (ip)  (p)-Bg(i/>)  (p)  |  <  bf  (p,q)-bg(p,q)  |<//(q)|dr~ 


<  1 1  *  I  L 1 1  £-9 1 1 1  1  =2j  IHI*war9 

ro 

This  establishes  the  Holder  continuity  of  f^B^.  It  then 
follows  from  the  assumed  continuity  of  the  map  f»-gf  that 
the  map  fn-Gf  from  Uad  into  C(rQ),  where  Gf  is  defined  by 

(2.19)  Gf (p) :»Bf (gf) (p)  =  j  bf (p,q) gf (q) dr~  , 

ro 

is  continuous.  We  may  now  rewrite  the  integral  equation 
(2.14),  using  the  definition  (2.19),  as 

(2.20)  $(p)+  af (P,q) $ (q)dr~  =  Gf(p),  perQ  . 

:  i 

0 

This  equation  is  of  exactly  the  same  form  as  that  treated 
in  Theorem  2.2  of  [  1  ].  The  remainder  of  the  proof 
follows,  almost  verbatim,  the  proof  in  [  1  ]  and  we  do  not 
repeat  the  arguments  here. 

The  continuity  of  the  mappings  f^$(*;f)  and  fH-G^ 
lead  immediately  to  the  following  corollary  which  establishes 
the  existence  of  solutions  of  our  optimization  problem. 


Corollary  2.1:  Under  the  hypotheses  of  Theorem  2.1,  the 

functional  J [ • ]  defined  by  the  equation  (2.8)  is  continuous 

as  a  map  from  Uad  into  S  and  consequently  takes  on  its 

absolute  minimum  on  the  set  U  . . 

ad 

There  is,  of  course,  no  guarantee  that  the  solution 
of  our  minimization  problem  is  unique.  We  can,  however, 


formulate  a  continuous  dependence  result  for  the  solution 
set  of  the  minimization  problem.  To  this  end,  we  consider 
the  functional  as  a  continuous  map  on  the  product  space 
C(rQ) xc(rQ) •  For  example,  we  may  think  of  the  added  mass 
functional  defined  in  equation  (2.8)  as  a  functional 
depending  on  the  given  data  g  as  well  as  on  f. 

For  a  given  functional  J(f,g)  we  set 

(2.21)  i  =  inf { J(f ,g) | f eUad> 
and  define  the  set-valued  mapping 

(2.22)  $(g)  =  {feuad| J(f,g)=ig>  . 

Note  that,  by  Theorem  1.2,  for  all  geC(rQ),  $(g)^0.  In  fact 
if  geVcC ( Tq) ,  V  closed,  then  the  graph  of  this  set-valued 
mapping  is  closed,  as  is  shown  in  the  following  theorem. 

Theorem  2.2:  Let  l/cC(rQ)  be  closed  and  J:Uad*V->-R  be 
continuous .  Then 

(a)  the  real-valued  function  gw-i  is  a  continuous 
map  from  1/  into  R,  and 

(b)  the  set-valued  function  <J>,  defined  by  the 
relation  (2.21),  has  a  closed  graph. 

Proof :  Suppose  that  { gn }  and  {fn>  are  sequences  in  V  and 
Uad  respectively  such  that  gn**-g  in  C(rQ)  and  f  -*-f  in 
C1' 1(rQ) ,  and  fne$(9n) •  Since  U&d  is  compact  V  is  closed, 
f£Uad/  and  <i>(g)  is  defined.  We  wish  to  show,  first,  that 
i  -*i,  and  second,  that  f€$(g). 
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Now,  let  fge$(g).  Then,  the  continuity  of  J  implies 
that  i  ,g)  .  Moreover,  by  definition  of  i  , 


'n 


J(f0,gn)>i^  for  all  n.  Again  by  continuity  of  J, 


'n 


J(f0,g  W(fQ,g)=i  .  Hence  i  <J(f,g)=lim  i  <lim  J(fQ,g  ) 

y  y  n-Ko  yn 

Therefore  J(f,g)=i  =lim  i  and  fe$(g). 

9  gn 


Corollary  2.2:  The  optimal  solutions  of  the  minimization 
problem  for  the  functional  (2.8)  depend  continuously  on 
the  boundary  data  g  in  the  sense  of  Theorem  2.2. 


Optimality  Conditions 


We  turn  now  to  the  question  of  finding  necessary 

conditions  for  the  constrained  optimization  problem. 

Specifically,  we  will  treat  the  problem  of  minimizing  the 

functional  (2.8)  representing  added  mass,  under  the  additional 

assumption  that  the  subset  of  admissible  surfaces,  Uad,  is 

convex  as  well  as  compact.  Such  an  additional  requirement 

is  not  overly  restrictive  since,  as  we  remarked  in  section  2, 

the  class  of  all  admissible  surfaces  is  convex  in  c^'^(ro). 

If,  as  suggested  there,  we  take  Uad  to  be  the  class  of  all 
2 

admissible  C  -surfaces,  Uad  will  be  convex. 

As  is  well  known,  the  basic  necessary  condition  in 
this  context  takes  the  form  of  a  variational  inequality. 

For  convenience,  we  recall  here  some  basic  definitions. 

Definition  3.1:  Let  X  be  a  vector  space,  Y  a  normal  linear 
space  and  A  a  mapping  defined  on  a  subset  DcX  with  range 
in  Y.  Let  Xq£D  and  <J>eX.  If  the  limit 

(3.1)  A'  (xQ ;  $)  =  lim  i  [A  (xQ+e  <J>) -A  (xQ )  ] 

€+0 

exists,  then  A  is  said  to  have  a  Gateaux  differential 
A'(Xg;<i>)  at  Xq  in  the  direction  0.  If  A'(Xq;4>)  exists  for 
all  <j>eX,  then  A  is  said  to  be  Gateaux  differentiable  at  Xq  . 

If  A'(xq;4>)  exists  for  all  <J>£X  and  if  the  map  4>h-a  '  (Xq  ;  <j>) 
is  linear  and  continuous  then  this  map  is  called  the  Gateaux 
derivative  of  A  at  the  point  xQ .  In  this  case  we  will 
denote  this  map  by  A'(xa).  Thus  A'(xn)  is  a  linear 


continuous  map  from  X  to  Y. 


Remark ;  We  note  that  it  makes  sense  to  consider  the  limit 

(3.1)  only  if  Xg+e<f>eD  for  all  e  sufficiently  small. 

The  basic  theorem,  stated  here  in  a  form  slightly  weaker 
than  usual  (compare  e.g.  [4,  p.119])  is  the  following: 

Theorem  3.1:  Suppose  X  is  a  vector  space,  K=X  is  a  convex 
subset,  and  A:K-+R.  If  XgeK  is  such  that  A(Xg)<A(x)  for  all 
xeK  and  if  A'(Xg;x-Xg)  exists  for  all  xeK  then 

(3.2)  A'(x0;x-Xg)  >  0  for  all  xeK  . 

Proof :  Since  xQ  minimizes  A  on  K  and  since  K  is  convex, 

Xg+e (x-Xg) eK  for  all  ee[0,l]  and  all  xeK  and 

A(Xq+€ (x-Xq) )  >  A(Xg)  . 

Hence 

A(Xg+€ (X-Xg) ) -A(Xg) 

€  “  0 

for  all  ee(0,l]  and  all  xeK.  Since,  by  hypothesis, 

A'  (Xg;x-Xg)  exists  for  all  xeK, 

A(x  +e(x-x  )-A(x  ) 

(3.3)  A'(xn;x-xn)  =  lim.  - - - - - —  >  0  . 

e  +  0 

Remark:  The  slight  difference  in  the  form  of  the  statement 

of  Theorem  3.1  from  that  usually  found  in  the  literature 
(e.g.  [4]  loc.cit.)  lies  in  our  assumption  that  A'(Xg;h) 

exists  only  for  those  h=(x-Xg),  xeK.  This  is  crucial  for 
our  considerations  since  the  functional  we  treat  is  only 
defined  on  a  convex  set  (C^'^(fn)  functions  representing 


surfaces  in  the  Maz'ja  class)  and  cannot  be  interpreted 
outside  of  this  convex  set.  Thus  it  is  impossible  for  us 
to  consider  arbitrary  '  ^"-perturbations  of  the  form  x^+eh. 

The  inequality  (3.2)  is  the  basic  variational  inequality 
which  expresses  a  necessary  condition  for  the  problem  of 
minimizing  A  over  K.  We  remark  that,  if  A  is  convex,  then 
the  condition  (3.2)  is  also  sufficient  for  a  (unique) 
minimum.  We  will  not  consider  this  question  in  what  follows. 

The  application  of  this  result  to  our  concrete  situation 
(x=C1'1(r())  ,K=Uad)  will  be  possible  provided  we  can  show  that  the 
limit  in  (3.3)  exists  for  the  functional  J  defined  in 
(2.8).  If,  indeed,  this  limit  exists  for  all  fj^ad 
i.e.  if,  for  every  fleUad' 

Jlf.+eff-f,)  1-Jtfj] 

J'  lf1?  f-f x]  :-lim+  - - - ^ - — 

£-*■0 

exists  for  all  feUad,  then  solutions  of  the  corresponding 
variational  inequality  will  be  candidates  for  optimal 
solutions.  Notice  that,  in  light  of  the  existence  theorem 
of  section  2,  there  is  no  question  concerning  the  existence 
of  at  least  one  solution  of  the  variational  inequality 

(3.4)  J’ [f^;f-f^]  >  0,  for  all  feUad  , 

since  this  inequality  will  necessarily  be  satisfied  at  any 
optimal  solution  in  Uad.  What  is  needed  here  is,  first, 
the  calculation  of  J*[f^;f-f^]  and,  second,  the  description 
of  an  appropriate  approximation  scheme  for  the  solution  of 
the  resulting  variational  inequality  (3.4) . 

With  regard  to  the  differentiability  of  J  we  remark 
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that  a  number  of  authors,  see  e.g.  [10],  [3],  [9],  [24], 

and  [4]  have  discussed  similar  questions  with  regard  to 
what  are  generally  referred  to  as  "domain  functions". 

Many  such  discussions  (see  e.g.  [9],  [14])  consider  two 
dimensional  problems  and  rely  heavily  on  conformal  mapping 
arguments.  Here,  we  return  to  Hadamard's  technique  (see  [3]) 
to  establish  the  form  of  J*  .  This  derivation  depends  on  a 


technical  lemma  which  we  state  below  and  whose  proof,  which 

is  somewhat  involved,  we  will  relegate  an  appendix.  Recall 

that  the  functional  J  depends  on  the  surface  f  through  the 

solution  of  the  boundary  value  problem  (2.3).  Actually 

we  consider  a  special  class  of  boundary  value  problems, 

namely  those  for  which  the  given  values  of  ‘the  normal 

derivative  on  T  are  in  fact  the  values  of  the  normal 

derivative  of  a  potential  function  which  is  defined  throughout 
2 

fi  x[h,0].  That  is,  the  boundary  condition  (l.ld)  is 
replaced  by 


9u  _  3u 
"5n  '"Sn” 


on  r 


2  A  2  ^ 

where  V  u=0  in  R  *(-h,0).  Of  course  u  will  not  satisfy  the 
boundary  and  radiation  conditions  appropriate  to  the  problem. 
In  this  sense  u  plays  the  role  of  the  incident  field  in 
scattering  problems.  In  this  vein  we  define 


T 

u  -  u+u 


which  behaves  as  a  total  field  and  satisfies  the  condition 


(3.5) 


0  on  r  . 
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Lemma  3.1;  Let  u  be  harmonic  and  let  u^  and  u2  be  solutions 
of  the  boundary  value  problem  (2.3)  corresponding  to  the 
admissible  surfaces  F^rff^)  and  "2=r^2^  with  and 
denoting  the  interior  and  exterior  respectively  of  I\  , 
i=l,2.  Let  u^=u^+u,  i=l,2.  If  we  define  AJ  by 

9u.  9u~ 

(3.6)  AJ  :=  U;L  dri-  u2  dr2 

< 

i  yi 

1  2 

then 

AJ  =  u  dl^-  u  +  Vu^*Vu2dV-  Vul*VuldV 

2  rx  D"\(D“nD“)  D”\(D"nD") 

or 

(3.7)  AJ  =  j  [ j VG | 2+(Vu^-Vu2) ]dV-  [|Vu|2+(Vu^-VuJ) ] dV 

D2\(D"nD2)  D"\(D~nD~) 

The  proof  of  this  lemma  is  given  in  Appendix  B. 

We  may  now  use  this  form  for  AJ  in  order  to  compute 
J'  [f1;f-f1]  for  f^eU.^. 

Theorem  3.2:  Let  f^eU^  and  consider  the  functional 

(3.8)  J[f]  =  Re  <J>(p;f)g(p)drp,  feUad  ' 

r  (f) 

where  <p(p,f)  is  a  solution  of  the  boundary  integral  equation 

(1.9)  appropriate  to  the  boundary  F(f).  Then,  for  all  feUa(j, 
the  Gateaux  differential  J'(f^?f-f^)  exists  and  has  the  form 

(3.9)  J'(f1;f-f1)  -  Rej  (f (p) -f 1 (p) ) (f 1 (p) ) 2{ [Vu ( f (p) p+pQ) ] 2 

r 

0  +  [n><VuT(f  (p)  P+PQ)  ]  2}dT~  . 

Proof :  For  ease  of  notation  we  write  v=f-f^  and  let 
T  =r(f.+€v) ,  Oscsl.  Then 


(3.10)  i[J(£1+ev)-J(f1) ]  =  i  Re)  [ | 7u| 2+(Vuj-Vu^) ]dv 


•De\(Dl(lDt) 


II 

D~\ (D“nD“) 


[|Vu|2+(Vu^*Vu^) ] dV 


=  7  *elh-X2]  * 


Consider  the  first  of  the  integrals  in  (3.10) .  Then, 
iterating  the  integral  we  may  write,  using  spherical 
coordinates. 


(3.11)  I 


» 

Is*  _ 


Vu | 2+ (Vu^*  Vu^) dV 


DeX(DlnD£) 


■ill  ■ 

A1i-R(9,4>) 


[|Vu|2+(Vu^«7u^) ] r2dr (do 


where  A1  is  a  measurable  subset  of  rQ,  R ( 0 , 4>)  is  an 
interval  of  r-values  and  do=sin9d9d<)> . 


(©,<!>) 


ev  ( 9  ,  <j>) 


D”\ (D^nD~) 


_ 

€  De\(D1nD£) 


Figure  3 


Using  the  mean  value  theorem  equation  (3.11)  becomes 
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*1  = 


f  f  ff1(e#<M+cv(e,<i>)  2 


AiLfi(6/4)) 


[  | Vu| 2+(Vu^*Vu^) ]r2dr|da 


=  e 


v  (0 ,  <|>)  [f1(0f<O)+X  (e,e,*)v(6,<fr)  ]2  [|Vu|2+(Vu*-Vu^)  ] 


do 


r=f ^+Xh 


where  0<X£ev.  Similarly,  we  arrive  at  the  same  expression 

for  I2  recognizing,  in  that  case,  that  f  ^  (6  ,$) +ev(0  ,<J>)  sf  ^  (9  ,<J>)  . 

Note  that  AiuA2=ro'  that  X=0{e)  and  that,  according  to 

Theorem  2.1,  the  solutions  of  the  integral  equations  depend 

continuously  in  the  C1'1  norm  on  the  functions  f.  Hence, 

T  T 

decomposing  Vu^  and  Vu£  into  normal  and  tangential 

3uT 

directions  and  recalling  that  =  0,  we  have 


[J(f1+€V)-J(f1)  ]  =  Re |  v(f1+ev)  2  [  (Vu)  2+Vu'r,Vu^'] 


dr( 

r=f ^+Xv 


Re 


vf2[(vG)2+(nxVuT) 2]drQ,  as  e-0 


which  completes  the  proof. 

Note  that  the  differential  is  linear  and  continuous 
with  respect  to  v=f-f1  and  that  the  variational  inequality 
takes  the  concrete  form 


(3.12)  Re 


(f (p)-f1(p) ) (f1(p) ]2{ [Vu(f (p)p+p0) ] 2 
°  +[nxVuT(f (p)p+pQ) ]2}dr~  >  0 


Having  the  explicit  form  of  the  Gateaux  differential 
we  may  indicate  here,  very  briefly,  the  questions  which 
must  be  addressed  in  introducing  a  concrete  numerical 
procedure  for  the  approximation  of  solutions  of  the 
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variational  inequality  (3.12).  Systematic  discussions  of 
the  general  problem  may  be  found,  e.g.,  in  the  book  of 
Cea  [4]  or  in  the  review  paper  of  Oden  and  Kikuchi  {21] . 

We  repeat  that  here  there  is  no  question  of  the  existence 
of  a  solution  of  (3.12)  in  light  of  Theorem  2.1.  To  find 
an  approximate  solution  of  (3.12),  let  {Sh^o<h<l  a 
dense  sequence  of  finite  dimensional  subspaces  of  C^'^(ro) 
and  let  Uad=uadn®h*  Here  as  usual,  h  denotes  an 
appropriate  index,  normally  the  mesh  parameter.  The 
sequence  ^uacj^0<h<l  a  of  subsets  of  C1,1(F0) 

approximating  in  some  sense  the  constraint  set  Uad.  The 
approximation  of  (3.12)  will  then  involve  seeking  a 
function  flh*uad  such  that 

(3.13)  J’lflh:flh-fhl  *  0,  tor  all  . 


Note  that,  since  Uad  is  compact,  UadnSh  is  a  compact 
subset  of  a  finite  dimensional  space.  Again  we  may  appeal 
to  continuity  of  J  on  Uad  and  hence  of  the  restriction  of 
J  to  uadr|Sjl  to  establish  the  existence  of  solutions  of 
these  approximating  inequalities. 

Since  Uad  is  a  subset  of  the  finite-dimensional  space 
Sh,  one  may  express  f^h  and  f^  in  the  form 


N 


N 


(3.14)  f  *  l  ot .  4> .  ,  f  =  E  6 .  <J>.  , 
xn  i  i  n  i=1  i  i 


if  denotes  the  collection  of  basis  functions 

spanning  S^,  where  and  8^  are  components  of  f^  and  f 
with  respect  to  the  basis  functions,  respectively.  By 


3 


substituting  (3.14)  into  (3.13),  and  using  the  linearity 


of  J' (f ;v)  in  v,  we  then  reduce  the  variational  inequality 

rvN 


problem  to  an  optimization  problem  in  B  :  Find  {ct^eK^ 
such  that 


N 


N 


(3.15)  Z  (ct.-8.)J'(  Z  >  0,  for  all  {B,}eK 

A  -L  J  J 


i=l 


j  =  l 


where 


Kh:-j{Si)5*N|V." 

1  =  1 


Then  several  standard  algorithms  in  use  in  the  theory  of 
constrained  optimization  problems  may  be  directly 
applicable  to  our  study.  Four  major  methods  of  this  type 
are  available.  These  are  the  classical  method  of 
successive  approximations  (or  fixed  point  methods) , 
pointwise  relaxation  methods,  penalty  methods  and  Lagrange 
multiplier  methods.  However,  it  is  understood  that  the 
success  of  any  of  the  aforementioned  schemes,  and  the 
questions  of  convergence  of  the  solutions  of  the 
optimization  problems ( 3 . 15)  to  a  solution  of  the  variational 
inequality  (3.12)  will  depend  heavily  on  detailed  analysis 
of  how  the  approximate  constraint  sets  K^,  and  the  subspace 
are  chosen,  as  well  as  on  establishing  G&rding-type 
estimates  for  the  operator  J’ .  We  shall  pursue  these 
questions  in  a  separate  communication. 
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APPENDIX  A: 

In  Lemma  2.1  we  introduced  estimates  in  order  to  prove 
that  the  operators  Af  and  B^,  defined  by  equations  (2.14) 
and  (2.15),  are  continuous  on  C(Tq)  and  that  the  mappings 
fw-Af  and  f^-B^  from  Uad  into  B(C(ro))  are  continuous.  Here, 
we  outline  the  calculations  needed  to  establish  those 
estimates . 

We  begin  by  remarking  that,  since  the  set  Uad  is 
bounded  in  the  C^^-norm,  the  elements  of  Uad  satisfy  a 
uniform  Lipschitz  condition.  Hence,  writing  p^=f(p)p+pg 
and  qf=f  (q)  q+pQ,  where  we  affix  the  suffix  on  pf  and  qf  to 
accentuate  the  dependence  on  f, 

(ad  | pf-qf I  <  |f  (p)  I  (p-q)  l  +  |f (p)-f (q)  |q  <  M|p-q| 

where  M  is  independent  of  f  (recall  that  |f  (p)  | <b  for  perQ)  . 

Moreover,  we  can  obtain  a  lower  estimate  on  | p-q |  as 
follows  if  we  recall  that  for  all  feUad,  f(p)>a  pointwise 
on  rQ: 

|pf-qf  |  >  I  [f  (p)  (P-q)  ]-((f  (q)-f  (P))q]  I 
>  | I f (p) I |p-q|- I f (q) -f (P) I | 

from  which  we  obtain 

( A2)  |pf-qf !+| f (q)-f (p) |  >  |f(p)||p-q|  . 

But  since  f(p)=|f(p)|,  we  have 


|f  (q)-f  (P)  | 


If <q) I " I f (P) ! 


|qf-p0l-|Pf-Pol 


< 


Pf-qf 
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which,  when  combined  with  (A2)  and  the  pointwise  lower 
bound  on  f,  yields 

(A3)  ||p-q|  <  |pf-qf|  • 

Turning  now  to  an  estimate  for  the  kernels  a^,  we 
recall  (see  the  relation  (1.4))  that  the  normal  derivative 
of  y  can  be  written  in  the  form 

3  n(q-)  •  (Pf-qf) 

(A4)  an-  tv(pf,q()I  =  -  — —  -y  +  Q(Pf-qf) 

q  2ir|pf-qf[ 

where  Q  is  bounded.  The  unit  normal  to  the  surface  has 
the  form 


n  (q)  =  n(q)/Jf  (0qf<frq) 

where  (0  ,  <J>  )  are  the  spherical  angles  associated  with  the 

4  4 

point  qef(f)  and  n(q)  is  the  normal  vector  to  the  surface 
T (f)  whose  components  are  computed  in  the  usual  way  from 
the  parametric  equations  for  the  surface  in  spherical 
coordinates.  Explicitly,  if  q=f  (6q,4>q)  q+pQ  and 
q=(sin8  cos4>  ,  sine  sin*  ,  cose  )  then 

'■i 


n  (<j)  -  a?  x  ■  *  “f  Tn —  sin9  0  -f  *  ,  $+f  sin0  q 

"q  3^q  30q  q  q  3^q  q  q 


where 


Sq  ~  If"  and  “  (l/sin«q)  §|-  , 
^  q 


(9q'V  - 


f[f2sin26q+(||-) 


+  ( 


3f 


2.2 
sin  9  ] 
<3 


and 


Since  the  surface  is  described  by  a  function  feC  '  (rQ)  , 


the  usual  Lyapunov-type  estimates  yield 
( A5)  |n(qf) • (pf-qf) |  <ot|pf-qfj2 


hence,  with  (A4) , 


I  In"  Y(Pf'qf)  I  -  Tpf-qff  +  lQ(pf'qf)(  • 

q  r  r 

With  the  inequality  (A3)  and  the  boundedness  of  Q,  this 
last  estimate  becomes 


I liT  I  S  TT^T 

q  I p-q I 

and,  since  the  Jacobian  J^(q)  is  also  bounded  (see  the 
estimate  (2.13)),  we  have 


Q 

la f  (P»q)  I  35  l-jjj-  Y  (p£,qf )  (q)  |  <  7~T7  * 

q  I p-q I 

Thus  the  kernels  af  are  weakly  singular. 

It  is  much  easier  to  see  that  the  kernels  b^  are 
weakly  singular  and  we  omit  the  simple  estimates. 

In  order  to  establish  the  first  estimate  (2.17)  of 
Lemma  2.1  we  use  the  following  estimate  from  [  1  ;p.  52]: 


(A6) 


nf(qf)«(qf-pf)-ng(qg) 


(V?V 


<8 


~  2(1-6) 

p-q  I 


f-g 


5 

1,1' 


0< 


where  n^(q^)  is  the  unit  normal  to  the  surface  T(f)  at  the  point 
qfeF(f)  and  ng(qg)  is  similarly  defined  on  the  surface 
r  (g)  . 

Now,  using  (Al) ,  (A3) ,  (A5) ,  and  (A6)  we  have 


4s» 


nf(qf)-(gf-Pf)  _ 


|q£-pfl 


'Vpg 


n£(qf)-(qf-Pf)-n  tq  )  •  (q  -p  > . 

- -  ♦  w  W 


l^f-Pf I  3  kg-pg|3/ 


-'i-  ~  ,-1-26 

<  c | p-q | 


I  I  f“Si  i  i^+cjp-ql  2 

t 1 gq-pq I 3- I g£-Pf I 31 [ I g„-pq I 3- I q£-p£ I 31 { 


p-q  1 6 


or,  since  both  lqg~Pgl  and  jqf-pf| 
(A7)  Asolp-ql"1"26 I |f-g| l£>1+d| p-q 


satisfy  (Al) , 


-4 I P-q I3*36 


iqg~Pg 


3-|qf-pf I 


The  last  factor  in  the  second  term  may  be  estimated  as 
follows  where  we  write  Af =qf-pf=f (q) q-f (p) p ,  and  similarly 
for  g: 


g(q)q-g(p)p| 3-|f (q)q-g(p)p| 3 


|Ag|-|Af | 


|Ag|2+|Ag| |Af |+|Af | 


=  A*B  . 


Then 

(A8)  A  <  | Ag-Af |  <  |g(q) -f (q)  |  + |g(p) -f (p)  |  <  2 |  l f-g !  I  1#  X 

while  each  term  in  B  is  of  the  form  | f (p) p-f (q) q | <c | p-q | 

2 

and  so  B<3c|p-q|  . 


Combining  this  last  result  with  (A7)  then  yields 
A  £  fc1|p-qi”1“25+c2lp-oj"(1+<S)  ]  |  |  f-g  (  x 


which  establishes  the  estimate  (2.17). 

In  order  to  establish  the  last  estimate  (2.18)  we 

write 


(A9)  |bf(p,q)-bg(p,q) |  =  ! Y (pf , qf ) Jf (q) -Y (Pg.qg) Jg (q) ! 

s  |v(pf,qf)-y(pg,qg) | |Jf(q) !+|y(pg,qg) I |Jf (q)-Jg(q) 

Since,  according  to  the  estimate  (2.13),  the  Jacobian  is 
bounded  we  may  estimate  the  first  term  on  the  right  hand 
side  of  the  inequality  (A9) ,  by 

B1t-|Y{pf.qf)-Y(pg,qg)||Jf(q)|  s  6  |  |-  -  I 

since  R  is  bounded.  Then 


\Jf-J  I  s  I f-g I [f-sin-e+f^+f^sin^ei^+lgl | [f  “sin^+f ‘  +  f ‘sin20  p 


<J>  0  J  \  ^>  \  |  L*-  ° J 

-  [g2sin20+f 2+f 2sin20 ]  ** 


■<p  0 ' 


f-gl l1/1+p2 


r  -  2  .  2Q , ,2 . r2  .  2Q ,  ,  2  .  20 ,  2,  2  .  2n, 

[f  sm  0+f,+fQsin  0]-[g  sm  0+g,+gQsin  0] 

<p  W  (p  u 


[f  sin  0+f z+f zsin  0] [g  sinz0+gz+ggSinz0 ] 


s  (P;l+3p2)  |  |  f-g  |  \lfl  . 


j  p 

Hence  B-  <  -  | | f-g | | ,  .  and  this,  combined  with  the 

|p-q| 

estimate  (A10) ,  yields 


(All)  |bf (p,q) -bg(p,q) 


1 1  f-g  I 


6 

1,1 


from  which  Equation  (2.18)  follows. 
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APPENDIX  B: 


This  appendix  is  devoted  to  the  proof  of  Lemma  3.1. 
The  functional  J(r)=juy^  dT  where  u  is  the  solution  of 

r 

the  boundary  value  problem  (2.3)  corresponding  to  the 
particular  surface  is  often  referred  to  as  a  domain 
functional.  As  we  recorded  in  the  introduction,  the 
consideration  of  optimization  problems  involving  such 
functionals  goes  back  to  Hadamard  [11]  .  Other  references 
have  been  cited  in  the  introduction.  Here  we  return  to 
the  strategy  of  Hadamard  in  order  to  establish  our  result. 


Proof  of  Lemma  3.1: 


We  begin  by  recalling  some  notation:  u  is  an  harmonic 
2 

function  in  S  x(-h,0)  while,  for  i=l,2,  I\  is  a  bounded 
surface  of  class  and  u^  is  the  corresponding  solution 

of  the  boundary  value  problem  (2.3),  i.e., 


(Bl)  Au^  =0  in  D.  , 


O  U  . 

( B 2 )  +  ku^  =  0  on  y=0  , 


U  KA  • 

(B3)  =  0  on  y=-h  , 


(B4)  37T  "  “IS  °n  Fi  ' 


(B5)  3 IT  -  ikoui  =  0(r_J5)  * 


We  divide  our  discussion  into  two  cases: 


which  we  will  encounter  meaningful.  Writing  u^=u^+u, 
1=1,2,  and  recognizing  that  au^/SnsO  on  and  that  Au1=0 
exterior  to  V ^  and,  in  particular  in  we  have,  with 

the  divergence  theorem, 

|  Vu^VUjdV  =  V  •  (u^Vu^)  dV  =  u2  dr 


/■  -  Su  f  afi  .  3u. 

=  J  u  Tn  dr+j  u2  JE  dr+J  u  Tn ~  dr+J  u2  7n~  * 
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Using  the  boundary  condition  (B4) , 


(B6) 


7“l‘Vu2dv  -  \  S  it  dr-|  “2  Tr r  dr+ 


12 


3u. 

~  1  , 
u  air  + 


u 


3u^ 
2  3rT 


But 


r  3u, 

<B7>  J  (S  ur  -  ui  l>r- 


<s  inr  -  “i  lt>ar  * 0 


since  u  and  u^  are  both  harmonic  in  D12*  (Note  carefully 
that  we  always  take  the  normal  to  as  the  exterior  normal. 
This  accounts  for  the  choice  of  signs  in  (B7).)  We  may 
now  use  (B7)  to  replace  the  third  integral  on  the  right  hand 
side  of  (B6) .  Hence 

3u_ 


(B8) 


T  T 

Vu^*  Vu^dV  * 


f  ~  3u  f 
U  ^  dF-J  U 


'12 


2  3n 


dr 


“i  IS 


-  3ul  3u  3ul 

u  3n"  dr_J  U1  3H  dr+J  U2  Tn~  dr  • 


But 

3Uj 


3u 


cTn-  "Tn 


on  r 


r  3u.  3u2 

i'  i=1'2'  and  J  <u2  -  U1  7n~)dr  =  0  * 


Indeed  this  last  relation  follows  since  both  u^  and  u2  are 
harmonic  outside  r2  and  we  have,  therefore, 


(u 


2  3n 


3u,  3u- 

-  ul  3TT)dr 


t  3u^  3u2 

=  "j  (u2  3n  U1  3rT 


)  dC. 


+  lim 
R-*-“ 


f  3ul 

(u2  TFT  '  U1 

■►00  * 


3u~ 

^r)ds 


3D. 
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where  3DR  is  the  surface  of  the  cylinder  DR  of  radius  R. 
Noting  that  u^  and  U2  satisfy  the  free  surface  condition 
(B2) ,  the  integral  over  the  free  surface  is  simply 

t  3ui  3u2  f 

J  (u2  JrT  -  U1  SrT)dCf  =  J  [«2  (-Itu^-u^-ku^  ]dCf  =  0  . 

Cf  Cf 

Moreover,  the  radiation  condition  insures  that  the  second 
term  on  the  right  is  zero.  For  detailed  calculations,  see 
Kleinman  [15;pp.  7-9].  Therefore  the  relation  (B8) 
reduces  to 


(B9 ) 


T  t 

Vul*Vu2dV 


;12 


3u, 

u.  ^zr-  df- 


or,  rearranging  terms. 


(BIO)  AJ  * 

r  3ui 

1  U1  TrT  dr~j 

* 

U2 

3u2  dr  -  ( 

vr  dT  - 

VuJ-VUjdV+f  |Vu| 2dV  . 

■1 

1 

-» 

2 

C12 

D12 

This  is  exactly  the  conclusion  of  Lemma  3.1  for  Case  I 
since,  in  this  case,  D^MD^d’)  ~Di2  and  D~\(D“nD~)  =0. 


Case  II:  The  surfaces  and  f2  intersect. 


B 

#7777Z7777ZV7777777777777777777777777T77777? 


Figure  5 


Here,  we  use  the  device  of  Hadamard  to  extend  the  result 
of  Case  I  (see  e.g.  Bergman  and  Schiffer  [3  ]).  Thus,  we 
introduce  an  auxiliary  surface  Tg  containing  both  r^  and  r 2 
in  its  interior  and  apply  the  results  of  the  previous  case. 
With  self-explanatory  notations,  we  may  write 


f  3ul 

(B11)  J  U1  TIT  dr‘J  u 


la  -W 


dr  = 


J  vu^-vu^dv+l  S  |H  ar-J  S 


(B12) 


*2  TrT 


■ 

7-  u 


a  5n 


dr  = 


Vu2*VUgdV+  u  dr- 


Subtracting  equation  (B12)  from  (Bll)  we  obtain 

’  3U.  3U«  t  m  rp  f  _ 

(B13)  j  U1  75T  ar+  u2  Sir  dr  =  Vu^vujdv-  Kvu> 

n  n 


'-fill 


dr  . 


Notice  that  the  right  hand  side  is  apparently  dependent  on 
the  surface  ra  and  the  solution  ua  while  the  left  hand  side 
is  not.  To  see  that  this  dependence  is  illusory,  we 
examine  the  first  two  integral  terms  on  the  right  hand  side 
of  (B13) .  By  appropriate  addition  and  subtraction  we  have 


c  > 
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(B14)  Vu^*Vu^dV-j  Vu^Vu^dV  = 


la 


2a 


T  n..T» 


VuJ«  ( Vu^-Vu^)  dV+ 


L  D 


la 


°1 


Vul*Vu2dV 


a 


Vu^*  (Vu^-Vu^)  dV+ 


a  1' 


( Vuj • Vu^) dV 


LD 


2a 


la 


=j  VuJ’* (Vu^-Vu^)dV-j  VU2* (Vu^-Vu^)dV 


la 


2a 


L 


D2\(D1nD2) 


vu^»  7u2dv- j 


T  T 

VuJ-Vu^dV  . 


D1V^(D1r'D2) 


But 


|  Vu^- (Vu^-Vu^)dV-j  Vu|* (Vu^-Vu^)dV 


la 


2a 


t  T  T  3u. 
(ua"u2)TiT  dF" 


a 


T  T  9u2 
(ua“ui)7n~  dF 


3u£  3u2  3u^ 

since  =  0  on  and  =  0  on  r2.  Also  =  0  on  Ta 


and  thus,  after  rearranging  terms,  the  relation  (B14)  becomes 


(B15) 


\  K 


• VuTdV- 


VuT*  VuTdV 

fc  21 


la 


T  T  T 

T  3  ,  T  T.  ,  T  T\  3Ua  ,  T  9U2  T  9U1  1 

ua  ^n(uru2)"(ul"u2)^rT  +  U1  Tn  u2  TrT  \ 


dr 


fT  t  r 

Vu1*Vu2dV- 

d'Md'od") 


T  T 

Vu^'VUjdV  . 


D1\(D1no2) 


Now,  writing  u!®u^+u,  i=a,l,2,  we  have 


(B16) 


J  |*a  3n(ul“u2)“(ul_u2)3F"  +  U1  3n  u2  TrT  J  1 

3  f  -  3  ,  . 3u  ,  3,  x 

*  u  TO  ul-u2  ~  ul”u2  9n  +  ua  TO,uru2> 


‘  <U1'U2,TO“  + 


T  T 

T  9u2  T  9ull 
U1  TFT  "  u2  3n”  dr  * 


An  argument  completely  analogous  to  that  given  above  in 


Case  I  allows  us  to  conclude  that 


J  ua  5n(ul‘u2>-(uru 

r  L 


2>  ?*r  dr  “  0 


and  so,  as  can  be  easily  checked  by  writing  out  all  terms 
the  right  hand  side  of  (B16)  reduces  to 


r  3u0  3u. 

J  (U1  7n"  “  u2  TrT)dT  =  0  * 


again  because  u^  and  u2  both  satisfy  the  free  surface  and 
radiation  conditions. 

Combining  this  last  result  with  equation  (B13)  and 
(B15)  we  have 


3u.  r  3u~ 

(B17)  U1  TtT  dr'J  u2  TrT  dT  =  ,  _ 


m  m 

Vul’Vu2dV 


D2\(D1hd2) 


7ul'Vu2dV'ff  G  TO  dr-[  G  TO  dr 

D][\(D^nD2)  r2  rx 


which  is  the  desired  result 
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